INFINITENESS OF ZERO MODES FOR THE PAULI OPERATOR 
WITH SINGULAR MAGNETIC FIELD 



GRIGORI ROZENBLUM AND NIKOLAI SHIROKOV 

Abstract. We establish that the PauU operator describing a spin-1/2 two- 
dimensional quantum system with a singular magnetic field has, under certain 
conditions, an infinite-dimensional space of zero modes, possibly, both spin-up 
and spin-down, moreover there is a spectral gap separating the zero eigenvalue 
from the rest of the spectrum. In particular, infiniteness takes place if the 
field has infinite fiux, which settles this previously unknown case of Aharonov- 
Casher theorem. 
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1. Introduction 

The presence of zero modes, eigenfunctions with zero eigenvalues, is a typical fea- 
ture for two-dimensional spin 1/2 quantum systems involving magnetic fields. Such 
eigenvalues were first found by Landau (see [15]) for the Pauli operator with con- 
stant magnetic field, and the multiplicity turned out to be infinite. Later Aharonov 
and Casher [3] calculated the number of zero modes for a bounded compactly sup- 
ported magnetic field, and this number turned out to be finite and determined by 
the total flux of the field. The conditions on the magnetic field were gradually re- 
laxed, see [16, 5], until in [8] the case of measure- valued magnetic fields was settled 
and an Aharonov-Casher type formula was established for a magnetic field being 
a regular measure with finite total variation, thus producing a finite number of 
zero modes. On the other hand, a weak perturbation of the constant magnetic 
field leaves the space of zero modes infinite-dimensional ([14]). In the paper by 
Shigekawa [19] it was established that if the field is sufficiently locally regular and 
separated from zero at infinity (or tends at infinity to zero sufficiently slowly) then, 
again, the space of zero modes is infinite-dimensional. On the other hand, an ex- 
ample in [8] had shown that if the total variation of the field is not finite, there 
may be no zero modes at all, even if the total fiux of the field, defined as a condi- 
tionally convergent integral, is nonzero. Under some rather restrictive conditions, 
infiniteness of zero modes was established for a periodic magnetic field, see [6], [7]. 

A new type of magnetic fields was recently considered in relation to the study 
of zero modes by Geyler and Grishanov in [9]. They have studied a system of 
equal Aharonov-Bohm magnetic solenoids placed at the points of an infinite double- 
periodical lattice in the plane. Neither of the previous results apply for this case, the 
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field being very singular and the total flux being infinite. Nevertiieless, the authors 
of [9] proved that such field produces an infinite-dimensional zero energy subspace. 
Moreover, both spin-up and spin-down null subspaces arc infinite-dimensional. This 
property is proved to be stable when one adds a constant (positive) magnetic field, 
of arbitrary size for the spin-down component, and not too large for the spin- up 
component. Further on, in [10] this result was extended to certain perturbations of 
this periodic structure. 

Not so much is known about the rest of the spectrum of the Pauli operator. For 
the Landau operator, with constant magnetic field, the spectrum consists of Landau 
levels, eigenvalues with infinite multiplicity placed at the points of an arithmetical 
progression. Under a weak perturbation of the field, these eigenvalues, except the 
lowest one, may split, producing a cluster of the discrete spectrum around the 
Landau levels (see [19]). Thus zero remain to be an isolated point of the spectrum. 
On the other hand, a weak magnetic field without a background constant field 
leaves the whole positive semi-axis belonging to the spectrum, so no spectral gap 
arises. If the magnetic field grows unboundedly at infinity, the whole spectrum, 
except zero, is discrete (see, again [19]). Under rather restricting conditions the 
presence of the spectral gap was established in [6] , [7] for a periodic field. However 
in a more or less general case this question is still open. 

In the present paper we study the zero modes of the Pauli operator with a non- 
regular magnetic field with an infinite total flux. The typical example of the fields 
in question is a, probably infinite, discrete configuration of AB solenoids on the 
background of a more regular magnetic field. The Pauli operator for a strongly 
singular field is not essentially self-adjoint, there is an ongoing discussion on which 
self-adjoint extension of the Pauli operator in the presence of AB solenoids more 
adequately describes the real physical situation ~ see [1], [4], [22] and references 
therein. It turns out that depending on which approximation to AB field by more 
regular fields is chosen, with simultaneous adjustment of some other physical pa- 
rameters, different self-adjoint extensions can arise. Our main analysis deals with 
the so called maximal extension. Its advantage is its invariance with respect to 
singular gauge transformations reducing the AB fluxes. We handle also another ex- 
tension considered in the paper [8], also gauge invariant, but with different spectral 
properties. We discuss the relations of these two extensions in Sect. 2, as well as 
describe the connections of the study of zero modes with problems in the theory of 
analytical functions. 

In Sect. 3 we find rather general conditions for the infiniteness of zero modes 
and for zero being an isolated point in the spectrum, with a possibility to estimate 
the size of the spectral gap. We start by settling the long-standing hypothesis 
(see the discussion in [8]) that a field of constant direction with infinite total flux 
produces infinitely many zero modes. Further on, we show that this infiniteness 
is preserved under addition of a field with different direction, having a finite flux. 
This establishes Aharonov-Casher theorem for the case of an infinite total flux of the 
field. We pass then to the case when this 'wrong' component may have an infinite 
flux. Here, we suppose that the flux of the field through any disk of a fixed size is 
non-negative, at least far enough from the origin, moreover the fiux of the averaged 
field is infinite. This requirement, together with some additional local conditions, 
grants that the spin-down zero subspace is infinite-dimensional. If, additionally, 
the above local fltixes are separated from zero, then zero is an isolated point of 
the spectrum of the Pauli operator. For regular fields, this condition prevents the 
infiniteness of the spin-up zero modes, since spin inversion corresponds to changing 
the sign of the field. However, if the discrete component of the field is large enough, 
in other words, if sufficiently many Aharonov-Bohm solenoids are present, then it 
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turns out that the main condition can be satisfied for the spin-up component of the 
'maximal' operator as well, so there are infinitely many spin-up zero modes too. We 
also explain how the results change when we pass to the self-adjoint Pauli operator 
considered in [8]. Here the situation with both spin- up and spin-down zero modes 
does not appear. We conclude Sect. 3, by some examples, in particular the case of 
a periodic and quasi-periodic magnetic field fits into the general approach, and the 
results of [6] and [7] are substantially extended. 

Further on, we pass to the situation when the general results arc not sufficient, 
since the main condition of positivity of local flux may be violated. Supposing 
that the magnetic field in question is a perturbation of some initial field where a 
quadratic lower estimate for the potential is known, we establish such quadratic 
estimate for the potential of the perturbed field, thus ensuring the infinitcness of 
zero modes, but, probably, without the spectral gap. Among others, the constant 
one, AB-latticc, probably, on the background of a constant field, a periodic or quasi- 
periodic field with some mild local regularity may serve as the unperturbed field. 
Admissible perturbations arc rather general, in particular they allow existence of 
arbitrary large regions on the plane with field having 'wrong' direction. In the end 
we discuss some examples where the perturbation theorems can be applied. The 
estimates for the potential obtained on this way of reasoning, may be useful in the 
further study of the perturbation of the Pauli operator by an electric field. 

The starting point of our study was an attempt to understand the possibility of 
perturbing the results of [9], by means of changing the intensities and positions of 
the AB solenoids. We thank V. Geyler who attracted our attention to this kind 
of problems. Further on, when it turned out that much more general situations 
can be taken care of, the proof of the crucial theorem 3.2 appeared in the process 
of discussions with F. Nazarov. We highly appreciate also the discussions with 
B. Berndtsson on the spectral gaps and with L. Erdos about the definition of the 
Pauli operator. The second author (N.Sh.) was supported by the stipend from 
the Swedish Royal Academy of Sciences. Both authors thank the Mittag-LefHer 
Institute for hospitality when the work on the paper was in its most active phase. 



2. Definition of the operator 

We identify the real two-dimensional space with co-ordinates x = {xi,X2) 
with the complex plane C, setting z = xi + 1x2', as usual, d = dg = {di + 182) /2, 
d = dz = {d\ — i52)/2,5j = d^^, and the Lebesgue measure will be denoted by dx. 

Formally, the Pauli operator in ^2(1^^), with gyro-magnetic ratio 5 = 2, is defined 
as the square of the Dirac operator ^ 

V = a- (-iV + A-) = {ai{-idi + Ai) + a2{-id2 + A2)). 

Here ai , 02 are the Pauli matrices, ^'^—i^^ o)''^^^(^i O^)' ^^'^ 

real functions, components of the magnetic •potential A = (Ai, ^2)- So, "p = = 

-(a-(V + U))2. 

Introducing the notations Hj = —id^ + Aj, Q± = Hi ± ill2, we can represent 
the Pauli operator V as 

f r+ \ / o_Q+ 

\ V- J \ Q+Q- 
Formally, the operators Q± are adjoint to each other. 




Different sign conventions are used in the literature. We follow the sign choice made in [8], 
which is the opposite to the one made in, say, [5] . 
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The magnetic field is defined as B = curl A = diA2—d2Ai , and it is considered in 
the classical physics as the only actual physical reality, the potential being merely 

a mathematical fiction. This is the fact also in the quantum physics, provided 
the magnetic field (and therefore the potential) are not too singular. The latter 
statement means that if for two magnetic potentials Ai,A2 the equality B = 
curlAi = curlA2 holds in the proper distributional sense then the corresponding 
Pauli operators arc gauge equivalent; there exists a real function (jj such that the 
multiplication by exp(z0) transforms one of the corresponding Pauli operators into 
another. To assign a rigorous meaning to the above statement, one has to define 
the Pauli operator as a self-adjoint operator in the Hilbcrt space i2(M^), with a 
certain domain. The gauge transformation, the multiplication by exp(z0), should 
transform differential expression in the proper way as well as transform the domain 
of one operator to the domain of the other one. 

The standard definition of the Pauli operator by means of quadratic forms re- 
quires A G i2.ioc and is described, for example in [8] or [21]. Having in mind the 
representation (2.1) of the Pauli operator, the quadratic form 

PaM= I \a-{-iV + A)i,fdx=\\Q+i,+ f+\\Q.^lj.f,^=(^ ) (2.2) 

is introduced. If the magnetic potential A is sufficiently regular, say, A G ^4,100 
curl A e i2,ioc,divA e i2,ioc, the quadratic form (2.2) can, actually, be obtained 
from the expression {Vtp,ip), € by means of the integration by parts, thus 
justifying the use of (2.2) for more singular potentials. So, for A G //2,10c one can 
chose some domain for the form pA, where this form is closed, and accept the self- 
adjoint operator corresponding to this form as the Pauli operator. Unlike the case 
of the magnetic Schrodingcr operator, where all reasonable choices of the domain of 
the form turn out to be equivalent (see [20]), for the Pauli operator such equivalence 
is not established. In [8] the authors argue that the choice of the maximal domain 
consisting of all functions ^ L2 for which (2.2) is finite, is physically reasonable 
since this corresponds to the states with finite energy. 

For magnetic fields possessing local singularities, the ones which we arc going 
to study further on, this description is not satisfactory, as it was explained in 
[8], and another approach, based upon the scalar potential, was proposed. We 
will use the same way of defining the operator, with certain modifications. In 
what follows, the magnetic field will be represented by a Borel signed measure 
jjL having locally finite variation. We suppose, moreover, that the support of the 
discrete part of the measure does not have finite accumulation points. To the 
field /i we associate a scalar potential ^{x), a solution of the equation A\l/ = ji 
in the sense of distributions. The corresponding vector potential is defined as 
A = (^1,^2) — sgrad^f = {82'^ , — i9i^'), again in the sense of distributions. 

The quadratic form (2.2), under certain regularity conditions, can be transformed 

to 

p[V'] =p+[V'+]+P-[*-] =4 j |a(e-*^+)|2e2*(ix + 4 j |a(e*V-)|'e-2*dx. (2.3) 

For a field ji with singularities, is is the form (2.3) that is used for defining the 
Pauli operator. The decomposition of the measure fi, = /Xcont + Mdisc leads to a 
similar decomposition of the potential, ^ = ^disc + ^cont • We will use the potential 
^'cont constructed in [8]. This is a function satisfying the equation A^'cont = 
Mcont in the sense of distributions. It is established in [8] that such a function 
exists and possesses certain regularity properties, in particular, exp(±2^'cont ) G 

-^l.loc 
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If only the continuous part of the measure is present, the natural domain for the 
form (2.3) consists of all functions tp± for which (2.3) is finite. Although this defi- 
nition is rather implicit, this domain possesses an easily describable core: the space 
of functions tp± for which e~*'!/'+ and e^^- are smooth functions with compact 
support (see [8]). 

We assume next that only the discrete part of the measure fj, = /Udisc is present, 
Mdisc = 27r ^ ax5{x - A), a; G (2.4) 

The support A of /^disc will be supposed to be a discrete set, without finite accu- 
mulation points, moreover, uniformly discrete: 

dist(A, A\X) >r°,r° X), for any A G A. (2.5) 

Each component of the discrete measure is an Aharonov-Bohm (AB) solenoid (see 
[2]) with flux 2na\ and intensity a\. We consider the case of one solenoid first. The 
AB magnetic potential corresponding to one term in (2.4). B = 2na.S{x — x^); = 

{x^, X2), with intensity a can be chosen as A(a;) = (^—a ^'''^f'-' ,a ^^^f^ ^ , r = |a;— a;°|. 

The corresponding scalar potential is ^{x) = ^'(2) = alnjz — (here = 
X1+X2, from now on, the complex picture is more convenient.) Thus the expressions 
exp(±5') have singularities of the form \z — z^\^" at z^. 

Formally, the Pauli operator with AB field admits gauge transformations. For 
an integer m, we set <j){z) = cxp(— 2marg(z — z")). Then the multiplication by 
(f) transforms the AB Pauli operator with intensity a to the one with intensity 
a + m. The potential correspondingly transforms as \E' t— > ^'|2; — ^^l"*. Whether this 
transformation is a unitary equivalence of operators depends on how the self-adjoint 
operator corresponding to the form (2.3) is defined. 

To make our description of such self-adjoint operators more precise, we introduce 
the following notations. In what follows, the notations d, d have the meaning of 
derivatives in the sense of the space of distributions ^'(R^). For a closed set E, 
we denote by Oe, Qe the derivatives in the sense of ^'{E? \ E). 

We explain now the way of defining the Pauli operator, proposed in [8]. For 
a G [—1/2,1/2), one accepts as a domain of the form (2.3) the space of such 
functions ip G L2{p?) for which the derivatives d(e~^ipj^) and d(e^ip-) (thus taken 
in the sense of distributions in ^'(M^)) are functions, and the form is finite: 

p[il)] < 00. (2.6) 

With such domain, which we denote here by DEv{ct), the form (2.3) is closed and 
defines the self-adjoint operator which we denote by Vev = 'Psvioi)- For a outside 
the above interval, the operator is defined by means of the gauge transformation. 
For a given a, let a* be the unique number in the interval [—1/2, 1/2) such that 
a* — a = m is an integer. Then the Pauli operator Vev {01) is defined as VEv{cf) = 
exp(im aTg{z—z'^))'P Evict*) exp(— zm arg(z— z*^)). With this definition, the operator 
is automatically gauge invariant. However, for a ^ [—1/2, 1/2), the description of 
the domain does not agree with (2.6). In fact, if the distributional d{e~^4'+) is a 
function, for ^' = |z — z^j" , the gauge transformation leads to the expression 9((z — 
z°)'"e~*V+)) which is not necessary a function, it may contain the (5-distributions 
and its derivatives. This might be considered as a minor inconvenience, however 
it leads to the unnatural absence of invariance of the number of zero modes under 
the change of sign of the magnetic field, as can be seen from the version of the 
Aharonov-Casher theorem in [8] (or, more easily, from the non-symmetry of the 
main interval [—1/2, 1/2), chosen arbitrarily - see [18] for more details). 
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We consider, along with the above operator, an alternative one. For a given a 
we define Vmax = 'Pmax(Q!) as the operator corresponding to the quadratic form 

M = P+max[V'+] + P-max[V'-] 

= 4 J |a{,oj(e-*^/.+ )|2e2*da; + 4 J |a|,o}(e*V-)Pe-2*dx. (2.7) 

defined on such functions tp that the derivatives in (2.7) (thus understood in the 
sense of ^'(R^ \ {2°})) are functions and pmaxl'^] is finite. The operator T'max is 
again gauge invariant (see, again, [18] for corresponding calculations). 

Both constructions can be carried over to the case of a finite or infinite system 
of AB solenoids placed at the points of a discrete set A of the plane, as in (2.4). 
We say that the vector- function A(a;) ~ {Ai{x), ^2(2;)) is a vector potential for the 
magnetic field (2.4) if /Xdisc = curl A in the sense of ^'(R^). The function 4'disc 
satisfying the Poisson equation A^fdisc = Mdisc is the scalar potential. As above, we 
define the quadratic form p by the expression in (2.3). The gauge transformations 
enable changing all intensities ax by arbitrary integers. If mx, A € A, is a collection 
of integers then the gauge transformation changing the intensity at the point A by 
2Trmx consists in the multiplication by a function W{z). The function W{z), as 
proposed in [8], equals ^|^|| , L{z) ~ F{z)G{z) where F{z) is an analytical function 
having zeros of order mx at the points A where mx > 0, and G{z) has zeros of order 
— TOa at the points A where mx < 0. For any collection of ax, by adding proper mx, 
one obtains the rerfitced intensities G [—1/2, 1/2), used for defining the operator. 
As in the case of a single AB solenoid, the operator Vsviioix, ^ G A}) is defined 
as the one gauge equivalent to 'Pev{{<^*xt ^ G the latter determined by the 

quadratic form (2.3) on all functions ip for which this form is finite. Alternatively, 
the maximal operator 7^max({Q;A) A G A}) is defined for any set of (non-integer) ax 
by the quadratic form pA in (2.7) on the functions for which this latter form is 
finite. 

Finally, when both discrete and continuous components of the measure fi are 
present, the discrete one located at the points A of a discrete set A, the Pauli 
operators Vev ^-nd "Pmax are defined in a similar way, with only difference that the 
scalar potential 5* is now the sum of the potentials ^'disc and ^'cont corresponding 
to the discrete and continuous parts of the measure fx. We do not touch upon the 
question on which of these operators (if any) describes the actual physical picture, 
we keep however in mind that for a continuous measure as a field these operators 
coincide. 

In the general situation it is hard to describe the domain of these two operators 
explicitly. However, more can be said about the null subspace of these operators, in 
other words, about the zero modes. Note, first of all, that the quadratic form of pmax 
is an extension of the form pEV- Since both forms are non-negative, this implies 
that KeT{VEv) C KeT{Vmax)- This can also be seen from the direct description of 
the zero modes. If a function 'tp = (tp+jip-) lies in the null subspace of the operator 
Vev or Pmax than tjj must annule the corresponding quadratic forms pEV, Pmax, 
which means 

a(e-*V-F) = 0, 9(e*V-) = (2.8) 

for Vev and 

aA(e-*V+) = 0, aA(e*V-) = (2.9) 

for Pmax- Both (2.8) and (2.9) mean that the function = 6"*^+ must be 
analytical, /_ = e^ip- must be anti-analytical, but on different sets. For the 
operator Vev-, by (2.8), these functions must be entire functions of variables z,z 
respectively. On the other hand, for the operator Pmax, by (2.8) these functions 
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may have poles at the points A G A, but not too strong ones, so that still after the 
multiplication by exp(±^'), they get into L2- 

To make things more concrete we suppose that from the very beginning the gauge 
transformation is made, so that all intensities a\ are in the interval (0, 1) for "Pmax- 
In this case, the function e* behaves as |z — A|"^ near the point A G A. Therefore the 
condition dzf- = 0, together with /_ cxp(— ^E") € L2 leads to anti-analyticity of /_ 
at the points of A as well: near a point A S A the function must be summable 
with weight having a singularity of the form \z — A|~^"^, therefore the possible 
singularity of /_ is removable. On the other hand, for the spin-up component ip^, 
the function /+ has to be holomorphic outside A but near the points of A it must 
belong to L2 with weight |2: — A^"', which, due to ax G (0, 1), allows /+ to have a 
simple pole at A. This asymmetry can be reversed by changing the normalization 
of the discrete part of the measure: by mean of a gauge transformation we can 
decrease all intensities by 1 thus arriving at the measure fi' with negative discrete 
part having the intensity a'^ = (ax — 1) G (—1,0) at the point A G A. Then it is 
for the spin-up component that the null subspace is generated by entire functions, 
and for the spin-down one by meromorphic functions with simple poles. We can, 
moreover, take the first normalization when studying the spin-down component 
and the second one for the spin- up component, thus only entire functions will be 
involved. With this last agreement accepted, the derivatives involved in the forms 
do not depend in the space of distributions where they are considered, so we can 
painlessly omit the corresponding subscripts in our notations. 

To compare the null subspaces of two operators under consideration, we suppose 
first that all intensities ax lie in (0, ^). In this case for the spin-down component 
V- null subspaces coincide, being in both cases generated by entire anti-analytical 
functions. The null subspace of ^'+max may be larger than the null subspace of 
^+EV: since the latter subspace may involve meromorphic functions /+ in addition 
to entire functions for V^ev- If some of the intensities ax lie in [— i,0), both 
'P+max and P-max may have zero modes generated by meromorphic functions, so 
both null subspaces may turn out to be larger than the ones for Vev- See, again, 
[18] for more detailed comparison of these two self-adjoint extensions. 

Now we make some more observations about the part 4'disc of the potential, the 
one responsible for the discrete part of the measure. It follows from the uniform 
discreteness condition (2.5) that the discrete set A has a density not higher than 
that of a regular lattice, more exactly, 

N{R) = tt{A G A, A < i?} = 0{R^), r^oo. (2.10) 

Consider the sum 

*disc(-2) = aolog|^|+ + (2.11) 

if A = does not belong to the set A, the first term in (2.11) is omitted. The 
series converges uniformly on any compact set in C not containing the points in 
A, moreover the Laplace operator can be applied term-wise, so (2.11) produces the 
required potential. 

The particular case of a special interest is the one of a purely discrete measure, 
a regular lattice, with all intensities equal, 

A = {Xmim2} = wimi + W2m2, Q!A^j^2 = a G (0, 1), (2.12) 
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where wi, W2 are complex numbers with non-real LOxjuJi. In for this configuration, as 
it was noticed in [9] the potential is closely related to the Weierstrass cr-function 




(2.13) 



so that ^q{z) = alog |cr(z)| can serve as a potential for the magnetic field (2.4). It 
was established in [17] that, the potential Vto possesses a very special property: 

^o{z) = a{n{i^z^) + m\zf + p{z)), (2.14) 

where is a certain coefficient determined by the lattice, m = S being the 
area of the elementary cell, and p a A-periodic function, with proper logarithmic 
singularities at the points of the lattice. We do not care about the value of 
explicitly given in [17, 9]. Anyway, the first summand in (2.14) is a harmonic 
function; we subtract it and for a regular lattice we will consider 

^o{z) = a{m\z\'' + p{z)), (2.15) 

with p determined by (2.14). With such potential ^Po, for any entire function f{z) 
subject to 1/(2)1 < Cexp(7|zp) with7 < am, the function ■)/;_ (z) = exp(— a(m|2;|^+ 
p{z)))f{z) belongs to Z/2- This observation made in [9] proves infiniteness of zero 
modes for T'-max- Passing to the operator V^, we make the gauge transforma- 
tion reducing all fluxes to 27r(a — 1). For the reduced operator, the function 
^Q+{z) = {a — l)(rn|2p + p{z)) serves as a potential, and thus any entire function 
/+ of the variable z produces the L2 zero mode ijj^ = cxp((a — l)(m|zp + p{z)))f^ 
of the operator T'+max- As it is explained above, if a lies in the interval (0,^), 
infiniteness of zero modes holds also for V-ev, and the operator V+ev, as follows 
easily from the properties of 5*0, has no zero modes. In the case a € [— 1,0) the 
spin-up and spin-down components change their roles. In Sect. 4, we will show that 
such estimates for the potential are preserved under certain types kinds of the per- 
turbations of the regular AB lattice, thus, in particular, providing us with examples 
of magnetic fields possessing arbitrarily large regions with 'wrong' direction of the 
field, but, nevertheless, with infinitely many zero modes. 

So, under our normalization conditions, the study of zero modes is reduced to 
the study of existence of entire functions which, after being multiplied by a certain 
weight get into L2- This study can be done by means of explicit estimates for the 
potential 5", like in [3, 5, 16, 8, 9] or by indirect methods, cf. [19]. We are going to 
combine both approaches. 

3. Zero modes and the spectral gap. Methods of the theory of 

subharmonic functions 

In this section we establish the infiniteness of zero modes under rather general 
conditions. We start by proving this for the magnetic field with constant sign and 
infinite flux, and then relax; the positivity restriction in different ways. Further on 
we find conditions for the existence of the spectral gap. 

We suppose that the general conditions on the measure as formulated in the 
previous section arc fulfilled. 

Theorem 3.1. Let p be non-negative locally finite Borel measure on C, with the 
support of the discrete part not having finite accumulation points, and the normal- 
ization agreements of the previous Section be fulfilled. Then the operators V-ev > 
T- max have infinitely many zero modes. 

As it is explained in the previous section, it is sufficient to establish the following 
fact about (anti-)analytical functions, which is valid for any non-negative measure. 
Theorem 3.1 follows from it, with an obvious replacement of p by 2p. 
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Theorem 3.2. Let be a subharmonic function such that for the measure ji = , 

/x(C) = / dn{z) = oo. (3.1) 
Jc 

Then the spaces of entire analytical and anti- analytical functions f such that 

[ |/(^)|2e-*da; < 00 (3.2) 
Jc 

are infinite-dimensional. 

Proof. Of course, it suffices to establish just one part, say, about analytical func- 
tions, which we arc going to do, for convenience of references. The strategy of 
proving the theorem is the following. We chose a sequence of points z„, n = 1, . . . , 
in a special way. For any given n a collection of functions fk, k = 1, . . . ,n satisfy- 
ing (3.2) will be constructed in such way that fk{zi) =0, k < I, fk{zk) 7^ 0. Such 
system of functions is, obviously, linearly independent. 

Wc denote by D(z, R) the open disc centered at z with radius R. All the measures 
in the proof are supposed to be non-negative. 

We take zi = 0. Then we chose Ri so that /i(D(0, i?i)) > 46. Then for each 
A: > 1 we find R^ so that ii{VLk) > 20, = D(0, i?^) \ D(0, i?fc-i), ^^i = D(0,i?i). 
This can be done due to infinitcness of (3.1). Then we take any Zh, strictly inside 
f2fe, not in the support of /idisc , k > 2. 

Now wc fix n and find a sufficiently small positive S such that the disks D(zfc, 6) 
= 1, . . . , n lie strictly inside respective ilk- Then we fix measure < jJi supported 
in fii such that /io(f^i) — 26. Denote by 4*0 the logarithmic potential of measure 
/xo, ^o{z) = (27r)"^ / In |z- w|(i/io(w). The function behaves as *o(z) ~ v l-^l 
as \z\ — > 00. 

We set further = ^ — '^q, /^i = — /xq > 0, so = /xi and, by our 
construction, 

/xi(fifc) > 20,fc= l,...,n. (3.3) 

Further on we chose measures Vk < IJ-i, k — 1, . . . ,n, supported in respective fife, and 
such that Vki^k) = 20, and denote by Uk the logarithmic potential of the measure 
Vk, with the same asymptotic behavior C/fe ~ In \z\ for large \z\. For a positive 
h < 6/2, we denote by z^j^. A; = 1, . . . , n, the measures supported in the respective 
disks T){zk, h) and coinciding there with times the Lebesgue measures, so that 
the logarithmic potentials U'^ of these measures have the same asymptotic behavior 
for large \z\ as Uk, Ul.{z) ^ In \z\. We denote by U^{z) the subharmonic function 



U\z) = ^,{z) + J2iULiz)-Uk{z)). (3.4) 

k=l 

By our choice of measures, Uk and Uj^ differ controllably for large z. In fact, 

U',iz)-Ukiz) = ^Jln\l-^\ {dul - dvk), (3.5) 

and since for \z\ > 2(i?„ + 5), \'w\ < R„. we have ln|l - w/z\ < ln2, ho \Ul.{z) - 
Uk{z)\ < 201n2(27r)~^ < 15. Adding up such estimates for all k, we obtain 

|C/''(z) - 1'i(z)| < 15n (3.6) 

for large \z\, \z\ > 2{Rn + S). Therefore, for any non-negative function v, and any 
R>2{Rn + S), 



L 



?;exp(-*i)da; < e^^" / w exp(-[/'*)da;, (3.7) 

R<\z\<1R Jr<\z\<2R 
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For z small, \z\ < 2{Rn + 6) but lying outside the disks D{zk,S), we note that 
each function Ul{z), being the logarithmic potential of a measure supported in 
the disk D(zfc,ft,), h < S/2, is bounded by some constant depending on 5 and 
Rn, \Ul.{z)\ < C{d,Rn)- The potential Uk, being the logarithmic potential of a 
measure supported in the disk D(0, is not necessarily bounded from below 
but it is bounded from above, again, by some constant depending on 6 and Rn, 
U'kiz) < C{6, i?„) for \z\ < 2(i?„ + S). This gives us 

= Y,iU'kiz)-Uk{z)) > -2nC{6,Rn), \z\ < 2{Rn+6),z^ UD(zfc,(5). 

(3.8) 

Next wc fix a function cp G C^(C) such that cp vanishes in the disks D(zfe,(5), 
k < n, but (fi{z) = 1 in D(z„,J). To the function dip we apply the theorem by 
Hormander, see [13], Theorem 4.4.2, on the solutions of 9-equation in weighed 
spaces: we find a function g = Qh solving the equation Bgh = dip such that 

9h\^7-^-^^^dx< [ lat^pe-^'rfx. (3.9) 



+ 

We recall now that d(p{z) = in the disks D{zk,S). Therefore the estimate (3.9) 
gives 

-u'' r 
\2 ^ J„ ^ \^,„\2-u'-. 



\9hV7T-rT^dx < / |9^re-^ dx. (3.10) 

(1 + kl) MuD(2,,5) 

Now we use the estimates (3.6) and (3.8) which enable us to replace in the right- 
hand side of (3.10) the weight exp(— [/'*) by exp(— \l/i). We obtain therefore the 
inequality 

/ lg^'' nlTm2 ^^ < ^ / |a<Ppe-*Mx = K (3.11) 

Jc (1 + kr) Jc\uD{z^,5) 

with a constant C depending on n, S, Rn but not depending on h. The left-hand 
side in (3.11) can be estimated from below for large i?, using (3.7), which gives 

J \gh\^e-'''dx<{l+R^fK, (3.12) 



/ 



9h W, '\ , dx<2K, (3.13) 



R<\z\<2R 

as well as 

e 

(1 + 

\z\>2{rn+5) 

SO that the weighted norms of gh over any annulus R < \z\ < 2R are bounded 
uniformly in h (of course, the bound may depend on R) . Recalling now that is 
a smooth function with compact support, we deduce from (3.12) that the weighted 
L2 norms oi gh — f over the annuli are bounded uniformly in h as well. 

We set fh = gh — ^1 dfh = dgh — dip = 0, thus fh is an entire function, moreover, 

1/2 / X 1/2 



/ \h\dx < J l/Hfe-'^dx J 

\z\<2R \R<|z|<2i? / \fl<|z|<; 



e*i(ix I <C{R)K 

R<\z\<2R \R<|z|<2i? / \fl<|z|<2fl 

(3.14) 

It follows from (3.14) that in any annulus R < \z\ < 2R the family of entire functions 
{fh}, < h < 6, has bounded Li-norms, therefore it has bounded C^-norms of any 
order in a smaller annulus (see Theorem 1.2.4 in [13]) and thus, by maximum 
principle, bounded C^-norms in any disk \z\ < R. Therefore, by Montel's theorem 
(see, e.g., [12], Theorem 15.2.5,) this family is compact with respect to uniform 
convergence on compacts: there exists a sequence hi ^ Q and an entire function 
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/ such that fhi converges to / uniformly on any compact. This implies that the 
sequence of functions ghi converges to g = f — uniformly on any compact. But 
now note that for a fixed k, the potential Ul{z) equals 

Ul{z) = -\nh + b{\z-zk\^-h')x^ 

TT 

inside the disk D(zfe, h), while all the other terms in , sec (3.4), do not depend on 
h or are uniformly bounded in h. Therefore exp(— J7'*(z)) has the order h~^^/'^ in 
D(za;, h), and since the constant in (3.11) is independent of h, the sequence ghi (zk) 
may only have as its limit value. So, g{zk) ~ and f{zk) = g{zk) + f{zk) equals 
for k < n and 1 for k = n. This function / is the one we arc looking for, because 
Jr \.f\^e^'^dx = Jj, |/pe~*i~*°da;, and the finiteness of the latter integral follows 
from the estimate 

g-*o(«) < (-gxp ^-^y" A*orfa;ln|^|^ < Cexp(-13/7rln |.2|) < C|^|-^ 

for large \z\, so that /j, l/pe-^dx < C J \f\^e-'^^{l + \z\'^)-'^dx, which is finite 
due to (3.13). So we have found the function /„. The functions fk, k < n are 
constructed in the same way, just the function (p has to be chosen to be equal 1 in 
the disk ^{zk, 6) and vanishing in ^{zk',6), k' ^ k. □ 

Having established Theorems 3.2 and 3.1 for a non-negative measure, we have as 
our next goal extending the results to measures having a negative part. The general 
requirement here is that the negative part ^_ is in a certain, each time concretely 
defined, sense weaker than the positive part For the rest of the section we 
suppose that the measure /i+ has infinite flux, /j,+ (M^) = /i_|-(C) = oo. 

Let us first discuss what may be the obstacle for an entire function / with finite 
/ \ f\^e~^^+dx to be quadratic summable with the weight e~^*++^*-, where ^- 
is a potential for It may turn out that grows at infinity, so a certain extra 
decay of / is required. On the other hand, the local singularities of the potential \E'_ 
can only be negative, and they would not cause any trouble since the introduction 
of the weight e^*" can only improve the convergence of the integral of 

So, the easiest result in this direction concerns the case when we can explicitly 
estimate the growth of ^+ and then take care of the corresponding term in the 
weight. 

Corollary 3.3. Suppose that fj,- has compact support. Then the statements of 
Theorems 3.2 and 3.1 hold for fj, = iJ,+ — i^-. 

Proof. Let ^-(C) 27r<I>,<l> > 0. Then the logarithmic potential ^'-(z) of ^_ 
grows at infinity as $ln |z|. Let TV be some integer larger than $. Take TV points 
Zi,Z2, ■ ■ . ,Zn such that each Zk is not a common zero for the space of entire 
functions / with finite J \f\'^e^'^^+dx. The latter can, surely, be achieved, and this, 
in particular, means that e^^*+ belongs to Li near Zk. The conditions f{zk) = 
0, k = 1, . . . ,N define a subspace =^jv of co-dimension N in =$f, so jSfjv is infinite- 
dimensional. Now fix a polynomial p(z) having simple zeros at the points Zk. 
The polynomial grows as |z|^ at infinity, therefore all functions of the form g = 
p{z)~^f, f G =$fjv are entire and have finite integral J \g\^e~^^+~^'^^- dx. □ 

Relaxing the condition of the compactness of supp/i_, we suppose only that 
/U_ (C) = 27r$ is finite. 

Corollary 3.4. Suppose that ^-{C) < oo. Then the statements of of Theorems 3.2 
and 3.1 hold for /x = /i+ — /x_ . 
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Proof. Wc choose the potential for the measure /x_ in the form: 

*-W = 7^/ ln\z-w\dn-{w)+^ [ lnfc^d/._H =C/o(z)+;7i(z), 

^'^ Jd(0,5) ^'^ Jc\T>{0,5) M 

(3.15) 

For l^l > 5 we spht Ui into two terms, Ui = U' + U" , where the first term 
corresponds to integration over the disk |w — 2;| < \z\/4: and the second one to 
the integration over the rest of the plane. In the disk the expression ^^^~|"' is 
smaller than 1, the integrand in (3.15) is negative, and U' < 0. To estimate the 
second term we note that € [2/3,6] for \w\ > 2\z\, and G [^jTyi f] 

for 5 < < 2\z\. Thus the integrand in U" is majorated by Ci + C2ln|^;| and 
therefore |[/"(z)| < (Ci + C2 In |z|)/i({|z| > 5}). A similar logarithmic estimate, 
with coefficient /x(D(0, 5)), holds for Uq for large \z\. Thus, as a whole, we have 
'^-{z) < Clog \z\, and the proof is concluded exactly as the one for Corollary 3.3. 

□ 

If the negative part of the measure fi is infinite, infiniteness of zero modes can still 
be established supposing that /i becomes non-negative after an averaging, however 
we need some additional local regularity conditions. 

Further on consider the following conditions for the signed measure /i = /U+ — 
M± > 0. 

Condition 3.5. There exist constants ro > and Oq G (0, 1) such that ijl+{T){z, ro)) < 
2-k6q for any disk inM? = €- with radius ro . 

Note that Condition 3.5 implies that if AB solenoids are present, their intensities 

lie in the interval (0, 6*0) for Vmax- 

In particular cases we also suppose that /i+ and/or /i_ satisfy 

Condition 3.6. There is a constant Ai and a radius Ri such that for any disk 
T>{z,Ri) 

J \ln \z - w\\diJ,±{w) < Ai. (3.16) 

D(z,fli) 

In particular. Condition 3.6 is satisfied if the measure /i± is absolutely continuous 
with respect to the Lebesgue measure dx and the corresponding densities belongs 
uniformly to Z/p,ioc for some p > 1. Note also that if this condition is fulfilled for 
some Ri it holds for any other with a different constant Ai. 

For a measure jj, satisfying Condition 3.5, we consider a potential ^'(z), a solution 
of the equation = jj,, as well as the potentials of the measures = 
IJ.±, ^ = 4"+ — The potential '5 (as well as is determined not uniquely 
but up to an arbitrary harmonic function. 

The first elementary fact we establish concerns measures satisfying Condition 3.6. 
Let X be a smooth non-negative function x G C^{T){0, R)),x = x{\^\) for some R, 
such that /q^q X = 1) a-^id we set ^'_r = \E' * X) so that A'ifn = = jj, * x- We 
set also *±,ij = *± * X) A*±,_r = /x±,fl = ij.±*x- 

Lemma 3.7. Suppose that the measure fj,- or iJ.+ satisfies Condition 3.6. Then 
there is a constant C = C{R,Ri,Ai) such that, with the corresponding sign ±, 

\^±{z)-^U^)\<C. (3.17) 

Proof. Denote by D the disk D(z, i?) and by D' the concentric disk with twice as 
large radius. Split the measure /x+ into the sum of the measure supported in D' 
and /t+ supported outside this disk. Correspondingly, the potential ^+ splits into 
the sum of = |J.'_^_*Go and = Go{z) = (27r)~^ In \z\. The function 
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^+ is harmonic in D', therefore ^+ * X = ^+ in the disk D, in particular, at the 
point z. The potentials ^'^ and * x are bounded, by (3.16). This proves the 
Lemma. □ 

Next we will study the potentials of measures for which satisfies Condition 3.5, 
and another component, /i_ satisfies Condition 3.6. Fix some R. For a fixed point 
zq denote by Dq, D, Di, D2 the disks with center at zq and radii, respectively, 
ro /2, To, To + R,rQ + 2R. We fix a non-negative moUifier x G Cq^(D{0, R)) as above 
and set '^r{z) = x * *±,fl = X * *±- 

We prove now our main local estimate. 

Proposition 3.8. For a fixed z, suppose that, Condition 3.5 is satisfied for and 

Condition 3.6 is satisfied for /i_ in AR-neighhorhood of z. Then there exist constants 
Co = Co{ro,R,9o,Ai), Ci = Ci{ro, R,eo, Ai), andC2 = C2{ro, R,Oo, Ai) such that 

J e-^''^'^\f{z)fdx<Co J e-2*«(^)|/(^)|2dx + Ciy e-^''''^'^\df{z)fdx (3.18) 

Do Di Di 

and 

J e2*-.«W|/(^)|2da; < <^2 y e'''^-^'^\f{z)\''dx (3.19) 

Di Di 

for any function f € L2(Di), as soon as the inequalities make sense. The derivative 
d in (3.18) can be replaced by d. 

Proof. For brevity, we prove the inequalities for the disks centered at the origin, 
noticing that the constants in all estimates below depend only on ro,R,6Q,A\. 
First, due to Lemma 3.7, we can restrict ourselves to a non-negative measure /i 
since the negative part of ji contributes to the estimates only with a constant factor 
when passing from e~* to e~*". 

We split the measure into the sum jj. = fj,' + v, so that fi' is supported in 
the disk D2 and v is supported outside this disk. The potential ^' splits into two 
parts, = 5*' + iJ , where 'i' = G * fi' is the Newton potential of the measure , 
G{z) = l/(27r) In \z/R\ , and H{z) = \l/(z) — '^'{z) is a harmonic function inside the 
disk D2. 

Correspondingly, the smoothcncd potential r splits into two terms, r = 
"^'j^ + Hr, where ^'^ = \E'' * X) Hr = H * x- Note that since H is harmonic inside 
D2, the functions H and Hr coincide inside Di. The function ^f'^ = /z' * G * x is 
bounded in Di, j^f^l < ci = ci{ro,R, Ai,6o). 

Let g{z) be a function, anti-analytical in Di, such that H{z) = —\n{\g{z)\). 
Then we have for any / 



e-2*W|/(.)|^dx= / e-2*'W|/(z)<?(z)|^dx, 
[ e-2*«(-)|/(^)|2da;= / e-^'''n(^^\f{z)g{z)fdx, 
f e-2*«(^)|5/(^)|2d:c= / e-2*«(^)|5(/(z)g(z))|2d:c. 



/Di 

So, denoting u = fg and taking into account the boundedness of ^'^(2;), we see 
that it is sufficient to establish the estimate 



e \u\''dx <C j (|m|2 + \du\'')dx (3.20) 

Do Di 

To prove (3.20), we split /x' into further two parts, /u' = /iq + Mi where /iq is 
supported in D and ^,l in Di \ D. The function = /ii * G is bounded in Dq. In 
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fact, the distance between points in Dq and in the support of 122 hes between ro/2 
and R + 2ro, therefore 

1*1(0)1 <max(| ln(ro/(2i?)|, I ln(i? + 2ro)/J?|)|M|(Di) <C. 

If /io = 0, the required inequaHty is now obvious. Otherwise, in order to estimate 
the contribution of ^'o = /Xq * G, we apply Jensen's inequality: 

2i,iI3)J^G{z-w)^ ^ f -2MD)G(.-w) M^) ^ p2»(D) / i i-»(D)/^ ^/^H 

- Jd m(D) ' ' /"(D) ■ 

So, for the left-hand side in (3.20) we have 

/ e-2*°(-)|n(z)|2dx<C / [/ |u(z)|2|z_u;|-''(D)/-dJ (3.21) 

In the inner integral we apply the Holder inequality, taking into account that 
/i(D)/7r < 2610 < 2: 

J Do 

(3.22) 

provided q G (1, 00)) is chosen so that 6oq' < 1, therefore the norm of \z — w\~^^^^/'" 
in (3.22) is finite. The second integration in (3.21) gives then 



/ 



e-''">^^^\u{z)\'dx<C\\u\\l^^^^^y (3.23) 



Do 

Finally we apply the Sobolev type embedding theorem in the disk Dq: 

IHIL,(D0)<^^(IHIL(D0 + I|5'^IIL(D0) 

and recall (3.23) and Lemma 3.7. This proves (3.20), and therefore (3.18). The 
inequality (3.19) follows immediately from Lemma 3.7. Obvious changes establish 
Lemma for d. □ 

De-localizing (3.18), (3.19) leads to the following fundamental fact. 

Proposition 3.9. Suppose that for all z the measure 11+ satisfies Condition 3.5, 
and jjL- satisfies Condition 3.6, "^{z) is a potential for the measure /i and is the 
smoothened potential = \l/ * %• Then, with some constants Cq, C\,C2 depending 
only on ro,Ri,6o, Ai 

J e-^^^'^fiz)\''dx<Co J e-2*«(^)|/(0)|2da; + Ciy e-^^''^'^\dfiz)\''dx, (3.24) 



and 



J e-2*«(^)|a/(0)|2da; < C2 ^ e-^'^^'^\df{z)\''dx, (3.25) 



as soon as the inequalities (3.24), resp., (3.25), make sense. Again, d can he re- 
placed by B. 

Proof. We take a covering of C by the disks D with radius ro/2 such that the 
concentric disks with radius R + ro/2 form a covering with finite multiplicity >c. 
Then we write the estimate (3.18) for each disk D, and sum these inequalities. 
This leads us to (3.24). To prove (3.25), consider the usual splitting ^ = 5'+ — 
The function "if^ is subharmonic, therefore exp(— 2\E'+^i{) < exp(— 2\E'+), so it is 
suflacient to establish 

J e2*-'«(^)|/i(0)|2da; < C2 J e^'^-^'^h{z)\''dx, 
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where h{z) = exp(— The latter inequaUty follows immediately from its 
localized version (3.19). □ 

Proposition 3.10. Suppose that the measure /Lt_ satisfies Condition 3.6 Then 

j e-2*«(")|/i(z)|2dx <C2 j e-2*(^)|/i(2)|2dx, (3.26) 
c c 
as soon as (3.26), makes sense. 

Proof. To prove (3.26), note that the function is subharmonic, therefore 
exp(— 2^'+_ii) < exp(— 2^*+), so it is sufficient to establish 

J e2*-'«(^^|/(^)|2rfa; < C2 J e''^-^'^f{z)\''dx, 

c c 

where f{z) = exp{—^+)\h{z)\. The latter inequality follows immediately from 
Lemma 3.7. □ 

Now we can establish our next theorem on zero modes. We recall that for the 
operator Vmax we accept normalization of the discrete part of the measure /x so 
that all intensities of AB-solenoids lie in (0,1), while for the operator Vev these 
intensities lie in [— i), with excluded. 

Theorem 3.11. Suppose that the measure /U_ satisfies Condition 3.6, ji satisfies 
Condition 3.5, moreover, for a certain Ro > , A{z) = /i(D(z, i?o)) > for 
\z\ large enough and J^A{z)dx = 00. Then the spin-down components of the Pauli 
operators V- max; V-ev have an infinite- dimensional null subspace. If, additionally 
a stronger condition 

A{z) > Ao > (3.27) 

is satisfied for all \z\ large enough then the point zero is an isolated point in the 
spectrum ofV-EV, 'P-ma.x o-nd the spectral gap above zero is estimated from below 
by some constant depending on ro, i?o, -Ri, ^0, ^0, ^i- 

Proof. Note first that due to our normalization agreement, see Sect. 2, the null 
subspace in both EV and max cases is generated by entire functions. Therefore 
we suppress the corresponding subscript. Let ^' be a potential for fi. The null 
subspace of the operator V- consists of the functions of the form u e 1/2 (C) such 

that u = cxp(— ^l')/, with an entire (anti-analytical) function f{z). Take some 
R > Rq and a mollifier xo supported in the disk D(0, R — Rq). Set X = Xo * XRo' 
where xn.o is the characteristic function of D(0,i?o)- Under the conditions of the 
first part of the theorem, the potential ^'^(z) = * x is subharmonic outside a 
compact set. If / is an entire function, the second term on the right-hand side in 
(3.24) vanishes, so we obtain the estimate 



-2*(-)|/(z)|2d:, <Coj e-2*«(-)|/(z)|2d:c. 



(3.28) 



This inequality implies that if for some entire function / the function exp(— ^/j)/ 
belongs to L2 then exp(— ^f)/ also belongs to L2. However the space of the functions 
satisfying the former condition is infinite-dimensional by Corollary 3.3. Therefore 
the space of entire functions / with exp{—^)f e £2 is also infinite-dimensional. 

To prove the existence of the spectral gap, we note first of all that under the 
conditions of the second part of the theorem we can, by increasing Rq, have A{z) > 
Ao/2 for all z, so, that is strictly subharmonic. Prom Hormander's theorem 
(Lemma 4.4.1 in [13]) on solutions of the 9-equation in weighted spaces (we apply 
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this Lemma with d replacing d), it follows that for any function g € L2 such that 
exp{—'^ii)9 G L2, there exists a solution / of the equation df = g such that 

J e-2*«(")|/(^)l''^^ < c J e-2*«(^)|(;(2)|2da;. (3.29) 
c c 
Substituting this inequality into (3.24), and using (3.26), we obtain 

/ e-2*W|/(^)|2da;<c' / e-^^^'^\g{z)\''dx. (3.30) 
Jc Jc 

We set e~*^^)/(z) = u{z), e~'^^^^g{z) = v{z), and recalling that g = df and <5_ = 
e~*9e* we get the inequality 

II^^IlL <C||^^llL>^^ = Q-", (3-31) 
where v is an arbitrary function in L2 and u is a certain function in L2 satisfying 
Q-U = V. Let h{z) be the projection of u onto the subspace Af in L2 orthogonal to 
all solutions w of the equation Q-W = 0; we still have v = Q-h. The left-hand side 
of (3.31) can only decrease if wc replace there u by h. So, on the subspace JV which 
is the spectral subspace of the Pauh operator V- = Q*_Q-, corresponding to the 
nonzero spectrum, the norm of the function h is majorated by p-[/i], the value of 
the quadratic form of the operator V- on the function h. This exactly means that 
the nonzero spectrum of V- is separated from zero, and the width of the spectral 
gap is controlled by the constant C in (3.31). □ 

The restrictions imposed in Theorem 3.11 on the negative part of /j. can be 
relaxed. 

Corollary 3.12. Suppose that the measure ji satisfies conditions of the first part of 
Theorem 3.11 and v is a finite non-negative measure. Then for the measure fi — u 
the zero energy subspace is infinite-dimensional. 

The Corollary follows immediately from the main theorem and Corollary 3.4. 

If there are arbitrarily large regions where the field vanishes, we cannot use the 
second part of Theorem 3.11 to estabhsh the presence of the spectral gap, even for a 
non-negative measure since such regions violate the condition (3.27). The following 
statement shows that if there are such regions then, actually, zero cannot be an 
isolated point in the spectrum. We suspect that the result is known to specialists, 
but we could not find a reference, so we present a proof, for the sake of completeness. 

Proposition 3.13. Suppose that for any R > there exists a disk D with radius 
R where the measure fi is zero. Then zero is not an isolated point in the spectrum 

of the Pauli operators V- and Vj^ . 

Note that we do not suppose anything about the nature of zero as a point in the 
spectrum. 

Proof. To justify the statement, we find, for any e > 0, a function (j) orthogonal 
to the null subspace of the operator V- (or V+) such that p-[</'], resp., p+[0] is 
smaller than e||(/)|p. To do this (for for example), we fix a non-trivial smooth 
function uq > 0, with compact support in the unit disk Di. For some constant C, 
the estimate ||wo||l2 1^ C'I|V"o||l2 holds. Now, let D be a disk with radius R > 
and center at zq, such that the restriction of /x to D is zero, and u be the function 
in Cq° obtained from uq by the dilation and shift, u{z) = uq{R~^{z — zq)). The 
function u satisfies 

\\u\\l,<CR-^\\Vu\\l, (3.32) 
Since the magnetic field is regular (in fact, it zero) on the support of u, the 
function u belongs to the domain of the operator V- 
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Denote by v the projection of u to the nuh space of the operator V-. The 
function v also belongs to the domain of the operator Si- (the function v is zero 

if, in particular, the null space is trivial.) So, (j) = u — v also belongs to the domain 
of V- and is orthogonal to Ker(P_); tj) is nontrivial since v, being a zero mode, 
cannot have compact support (unless v = 0). The norm of the function cj) in 7^2(1^^) 
is not greater than the norm of u in L2(IR^), or, what is the same, than ||m||i,2(d)- 
At the same time, for the quadratic form of the operator V- , we have 

P4<P]=P4u]= I exp(-2*)|5,(exp(*)w)|2da;. 

However the magnetic field /i = vanishes in D, so ^' is a harmonic function in 
D. Taking into account that u has its support in D, by means of the usual partial 
integration we obtain that p[u] — |Vtipdx. Now, the inequality (3.32) gives us 

\mUm^)<cR-'P-[<l>], (3.33) 

for the function (p orthogonal to the null subspace of V. Supposing that there are 
disks of arbitrary size R, not intersecting the support of fi we obtain that in the 
neighborhood of the point zero there are infinitely many points of the spectrum of 
V-. □ 

The reasoning in the proof of Theorem 3.11 does not apply directly to the oper- 
ator V+max, if, for example, the measure /x is continuous. In fact, when we pass to 
V+, we have to replace /i by — /i and d by d. The latter change is not that essential. 
However the measure —fj, does not satisfy the condition —ii(D{z, Rq)) > 0, this 
quantity is negative and the line of reasoning breaks down in several places. 

The case when the game can be saved for the operator Pmax is the one with 
a measure fi consisting of the system of Aharonov-Bohm solenoids placed at the 
points of a discrete set A, with intensities a\, a\ S (0, 1), on the background of a 
continuous field f such that both and V- satisfy Condition 3.6: 

M = MAS + 1^, IJ-AB = 2n^ axS{z - A). (3.34) 

In this case, as explained in Sect. 2, the measure —fj. can be reduced by a gauge 
transform to the one containing AB-fluxes with intensities (1 — a\), i.e., to the 
measure 

fi = p,AB-i', fiAB = 27ry^^(l - ax)5{z - A). (3.35) 

AeA 

If the set A is infinite, does not have large gaps, and the numbers ai are separated 
both from and 1, it is possible that both measures /x, /i satisfy conditions of 
Theorem 3.11. This leads to infiniteness of zero modes and the spectral gap for 

both P+ max and V- max 

We formulate a special case, where Conditions on the measures are expressed in 
more geometrical terms: 

Condition 3.14. There exist positive numbers r^ ,Rq such that any disc 'D{z,Rq) 

contains at least one point in A, any disk D(A,r''), A G A, contains no points in A, 
other than \, and all intensities ax, A € A, satisfy olx € (^o, 1 — ^o) for a certain 

eo>o. 

Corollary 3.15. Suppose that the measure /x has the form (3.34) its discrete part 

satisfies Condition 3.14- Suppose also that both positive and negative parts of the 
continuous part v in (3.34) satisfy Condition 3.6, moreover, for \z\ large enough, 
^0 > ^'±(D(^, i?o)). Then both operators Pi, max have an infinite- dimensional null 
subspace. If, moreover, 6o > i'±{'D{z, Rq)) > > 0, both operators possess a 
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spectral gap. The size of the gap is determined by the numbers r'^ , Ro,0o, Ai, and 
do- 

The Corollary above covers, among other purely discrete measure, i.e., 

an infinite configuration of AB solenoids satisfying Condition 3.14, as well as such a 
configuration on the background of a constant or 'almost constant' magnetic field. 
For a regular lattice A and equal intensities and a constant background field, the 
infiniteness of zero modes was established in [9], [10]. 

On the other hand, if the measure /x is continuous (thus "Pmax = 'Pbf), so that 
there are no AB solenoids, we can establish the spectral gap also for V+. 

Corollary 3.16. Let the measure ^ satisfy Condition 3.5, /i(D(z,i?o)) > > 
for \z\ large enough, and both fi± satisfy Condition 3.6. Then the operator V+ has 
no zero modes and possesses a spectral gap. 

Proof. Let ^ be the potential of the measure /it, ^'/{ its averaging, so that = 

II * X- is a measure absolutely continuous with respect to Lcbcsguc measure with 
positive density separated from zero. As it follows from Proposition 3.10 and the 
conditions imposed on n, for any function /, 

J exp(2*)|/|2rfx >C J cxp(2*fl)|/|2dx. (3.36) 

Now for the smooth function n we use the commutational relation which gives 

J |(e*«(")«)a(e-*«(")M)|2dx = y" e-2*«(^)|5j(e*«(")«)|2da;+4 J \u\''diiR, (3.37) 
c c c 

And this establishes our statement. □ 

An important special case of the above considerations concerns periodic magnetic 
fields. Such configurations attracted interest in early 80-s. In the papers [6] and 
[7] for the case of a rational flux of the fleld over an elementary cell of the lattice 
the infiniteness of zero modes was proved as well as the existence of the spectral 
gap, and nothing has been done since. We show below that the above restriction is 
irrelevant. 

Corollary 3.17. Let the measure n be periodic with respect to some lattice in the 
plane. Suppose that Condition 3. 6 is satisfied for both positive and negative parts 

of the measure /i and, moreover, the measure of one cell F of the lattice is positive. 
Then there are infinitely many zero modes for the Pauli operator V and a spectral 
gap. 

In fact, for a periodic measure, the Conditions 3.5 and positivity of /i(D(z, i?o)) 
for large \z\ obviously follow from the positivity of the measure of the cell, and this, 
by Theorem 3.11, establishes both infiniteness of zero modes and spectral gap for 
The presence of a spectral gap for V+ follows from Corollary 3.16. 

For a periodic magnetic field the infiniteness of zero modes may be proved also in 
the following way. Let A be the lattice of periods of the measure /z, with the elemen- 
tary cell F ■ Let a{z) be the Weicrstrass function of the lattice A defined in (2.13). 
As it is explained in Sect. 2, with proper y, the function ^I'o = ln|(T(z)| — ^{vz^) 
equals m\z\^ + p{z), m = with a periodic function p{z) having singularities of 
the form (27r)~^ In \z — A| near each point A of the lattice A, |F| being the area of 
the elementary cell. Consider the potential 

^'^(2;) = J exp{'ifo{z-w))diJ,{w) = j \z-w\^diJ,{w) + j p{z-w)diJ,{w), (3.38) 

with integration over the elementary cell F of A. Since A^o = '^S{z — A), the 
function is a potential for the measure /j.. Prom the Condition 3.6 for /U±, it 
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follows that the second term in (3.38) is a periodic bounded function and positivity 
of the flux $ over the elementary cell produces the growth of the first term in 
(3.38) as $m|zp. Therefore any function having the form u = exp(— ^'^(z))/ with 
an entire function f{z) growing not faster that exp(c|2;p), c < <i?m, is a zero mode 
for the periodic magnetic field u. 

More generally, we can consider a quasi-periodic magnetic field. Let /i = ^ /Xfc, 
k — 1, . . . ,N , so that each measure /x^ is periodic with respect to its own lattice 
Afc. Suppose also that Condition 3.6 is satisfied for each measure /i^, and the sum 
$ = Ylk frfcT' ^^^^ being the area of the cell of the lattice of periods for /ij., is 
positive. Then the reasoning used for a single periodic field goes through. Again, 
the smoothened field /x * x is a measure with a bounded positive density, which 
guarantees, as in Corollary 3.17, infinitcncss of zero modes and the presence of 
a spectral gap. Moreover, the explicitly constructed potential, being the sum of 
potentials grows as c\z\^. 

Of course if the sum is negative, then the potential for the measure /i = /^fe 
is majorated by — c|zp, and this leads to infinitely many zero modes for the Pauli 
operator V+. 

We note here that the fields above can be, by Corollary 3.12, perturbed by any 
finite measure thus preserving the infiniteness of zero modes. 

In the following sections we consider more strong perturbations of the field config- 
urations described above. These perturbation preserve the quadratic lower estimate 
for the potential, thus guaranteeing infiniteness of zero modes, however we will see 
that they may destroy the spectral gap. 

4. Perturbations of the field and zero modes 

From now on we consider a more special situation than in the previous section. 
Suppose that for a certain magnetic field jUq we know that there exists a potential 
"ifo satisfying the growth condition ^'0(2:) > 7|zp for \z\ large enough, with some 
7 > 0. Then any entire analytical function / growing at infinity slower than 
exp(7 — e)\z\'^ satisfies /exp(— ^f) G L2, and thus the null subspace of the Pauli 
operator V- is infinite-dimensional. Of course, this does not guarantee the presence 
of the spectral gap, and, moreover, the spectral gap may fail to exist, as we show 
below. Now, let /xi be another magnetic field with a potential ^'i [z) satisfying the 
estimate *o(z) - 'I'i(^) = o(|zp) or even \^a(z) - ^i{z)\ < 7*|zp with 7* < 7, for 
large \z\. Then, of course, the potential 4*1 grows at infinity sufiiciently fast so that 
the infiniteness of the number of zero modes for the Pauli operator with field /Ui is 
granted. More generally, the estimate for ^"0 and the inequality for the difference 
'^o{z) — ^'1(2;) may contain some singular terms, as it happens in the presence of 
A-B solenoids. 

There are several types of the magnetic fields that we can take as the unperturbed 
field jjLQ. One obvious example is the constant magnetic field described by the 
measure dfjbQ = Bdx, proportional to the Lebesgue measure - here the potential has 
the form ^'o = -jl^p. As shown in the end of the previous section, a more general 
field, a periodic field with mild local regularity properties (Condition 3.6) with 
nonzero flux through the cell, or the even more general one, a quasi-periodic field 
with nonzero j^^, also possess potentials subject to required growth conditions. 

Another kind of the starting point of our study can the observation made in [9] . 
Let A = Ao be a regular, periodic lattice with periods u)i,ui2 and all intensities 
are equal, a\ — a E (0,1), A G A. The scalar potential ^!o+{z) is defined in 
(2.15). Since a lies between and 1, the function u{z) = f{z)exp{—'^o+{z)) = 
f{z)exp{—am\z\^)exp{—ap{z)) belongs to L2 as soon as the entire function v{z) 
grows at infinity not faster than exp7'|^;|^, 7' < am. Of course, there are a lot 
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of such entire functions v{z), in particular, all polynomials fit. This proves the 
infiniteness of zero modes for "P+max or, under the condition a G (0, 5), for V+ev- 
Taking into account the possibility of applying the gauge transformations discussed 
in the end of Sect. 2, the same reasoning applies to the operator "P-max- Again we 
can perturb the (discrete now) measure ^0 by a measure having a potential with a 
slower growth than ^I^o and with controlled logarithmic singularities. 

So, in order to determine which perturbations of the initial field /xq preserve 
infiniteness of zero modes, one has to know which measures possess potentials with 
prescribed control over the behavior at infinity and, if needed, at singular points. 
A number of such results exist in the literature, sec, e.g., [11], however they are 
not sufficient in our situation since they require an extra decay of the measure at 
infinity, the condition we aim to avoid. Moreover, they do not usually take into 
account the possible cancellation of the contribution of the positive and negative 
parts of the perturbing measure. In this section we present some results on the 
estimates for potentials for certain classes of measures. 

The situation we consider first is the one when the whole (or a part of) mea- 
sure fiQ is re-arranged, more exactly, this measure is replaced by its image under 
some mapping of the plane. Under such re- arrangement, large regions with field 
having 'wrong direction' may arise, so that the positivity conditions of the general 
theorems are substantially broken. Nevertheless, the infiniteness of zero modes is 
preserved. We consider the case of /i being a continuous measure with fj,± satisfying 
Condition 3.6. Let : — * be a Borel measurable mapping and /x* be the 
measure induced by this mapping: for a Borel set E, ijl*{E) equals n{^~^{E)). 

Proposition 4.1. Suppose that the mapping $ satisfies the condition 

|$(w) - A| < a|u;r, r < 1 (4.1) 

for some a > 0, for \w\ large enough. Let also satisfy Condition 3.6. Fix some 
R> 0, such that (4.1) is satisfied for \w\ > R and define the function 

^*{z)= j z/^{w))-\n{l- z/w)]dn{w) + 

\w\<R 

j ^\}.n{l-z/^{w))-\n{l-z/w) + z{^{w)-'^ -w-'^)\dn{w) (4.2) 

|u;|>iJ 

Then the function ^* satisfies the equation A^* = 27r(/i* — and, moreover 

\^*{z)\<C\z\^+^. (4.3) 

Proof. The fact that the integral converges is obvious. To check that it satisfies the 
Poisson equation, it is sufficient to notice that for a bounded domain Q. containing 
the point z, 

j mn{l- z/^{w))dn{w) = j mn{l- z/w')dn*{w'). (4.4) 

To prove the crucial inequality (4.3) we have to estimate only the second integral 
in (4.2) since the first one may grow at most logarithmically at infinity. 
For a fixed z, we split the plane into three regions 

9.o = {weC, R< \w\ < |^;|/2},0i = e C, 1^1/2 < |«;| < 2\z\}, 

0.2 = {we C, 2\z\ < \w\}. 

and estimate the corresponding integrals Iq, /i, /2 separately. 



ZERO MODES OF THE PAULI OPERATOR 



21 



For the integrand in I2, wc have 

z 



1^ f 1 - ^) - In (1 - -) + ^ --]= (4-5) 



In 1 



In 1 



z 
w 



z 

2^ 



<f)(ui) 



and correspondingly, I2 sphts into I21 and l22- For \w\ > 2\z\, supposing |2| is large 
enough, we have |z/$(w)| < 3/4. 

Consider the function /i(^) — ln(l — + ^ + analytical in the open unit 

disk. For |^i|, |^2| < |, wc have 



- < la - 61 max < 4|6 - 61 maxdai', 161')- 

te[{i,{2j 

We set here ^1 = 6 = ^, obtaining from (4.6) 

z 



(4.6) 



h 



$(w) 



< 4 



$(w) 
1 



$(w) 



Therefore we get 



^21 



<C\z\ 



z 
w 

1 

w 

-4+T 



max 



$(w) 



w;r*+^d|M|(«^) < C\z\ 



1 + T 



(4.7) 



to|>2|z| 

Further on, in order to estimate I22, we make the transformation 



2 Vu;; 



which gives 



and therefore 



2 \w) 



1 
2 



1 2/ ,,w; + $(w) 



-3+T 



1/22(^)1 < C\z\'' [ \w\-^+^d\tx\{w) < Ca\z\ 

Taken together, (4.7) and (4.8) give 

1/2(^)1 < C\z\'+'. 
We pass to estimating Ii . We write the integrand now as 
SR [ln(l - z/^w)) - ln(l - z/w) + 0($(w)-^ - w' 



l+r 



In \^{w)/w\ + In 



$(w) 



(4.8) 
(4.9) 

(4.10) 



In the first term on the right in (4.10), we have |1 — ^{w)/w\ < 1/2, therefore 
|ln|$(w)/w|| = I In 1(1 - ($(w) - ■w)/w)\\ < C|w|~^+'^, and we get the estimate 
by for the integral over fii. The integral over fli of the second term in 

(4.10) we split into the sum of In and /12, the former being the integral over the 
domain fin : \w — z\ > 40121"^ and the latter over Q12 — iln. For In, we have 
< i, so 

liu— z| — 2' 

|(f>(w) - w\ 



In <C 



■/nil F - ^1 inn 
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w—z 



< |ln|$(u') 



For the integral over f7i2 we apply the rough estimate | In 

I In |w — 2;| I and integrate each term separately. Since the integration here is per- 
formed over the disk with radius 40121"^, the conditions on the measures fj, and fi* 

imply that I12 can be estimated by |zp'^|ln|z|| = 0(|z|-'-+'^). 



Finally, for Iq note that for to G f2o we have 
sufficiently large, therefore 



\Io{z)\< [ |ln|$H«;-i||d|/.|(«;)+ / |ln| 

Jqo Jno ^~ ^ 



w — z 



< 



2a|wr 



< 1/2 for \z\ 



+ \z\ [ \w-^ - ^w)-^\d\fi\{w) <C [ 
+C [ \w\^\w-z\-^d\fi\{w)+C\z\ [ \w\ 
All integrals in this expression are estimated by C^l^;]^"'"'^ 



w\-^+''d\fi\{w) 



-l+T 



d\n\{w). 



(4.11) 

□ 



Following the general pattern described in the beginning of the Section, we arrive 
at the following case where the infiniteness of zero modes is granted. Note that the 
somewhat complicated conditions imposed on the unperturbed measure are aimed 
to cover the interesting cases described in the begirming of the Section. 

Theorem 4.2. Let /xq be a measure, with discrete part Ho^disc = X^agAq '^a<5(-2 — A) 
supported on the set Aq satisfying (2.5), with Q < a\ < 6q < 1 for the operator 

'Pmax o^nd a\ € [—1/2, 1/2) for Vev ■ Let Ao(z) he the point in Ao closest to z (or 
any of such points) and oco{z) he ocx^(^z)- Suppose that the measure jJo admits a 
potential ^0 subject to the estimate 



^oiz)-a{z)\n\z-Xo{z)\ >7kl',7>0 



(4.12) 



for sufficiently large \z\ (if A is empty, In |z — A(z)| in (4.12) is replaced by zero). Let 
jjL be another measure and ^ be a Borel measurable mapping so that the conditions 
of Proposition 4-1 ire satisfied. Then for the measure hq + ji — n* there exists a 
potential ^* also satisfying (4.12) with some 7',0 < 7' < 7, and thus there are 
infinitely many zero modes for the Pauli operator V-max, resp., V-ev with the 
field /.iQ + 11 — . 

The version of Proposition 4.1 for the case of a discrete perturbing measure /i is 
also valid. We only give the formulation here, the proof being practically the same 
as above. 

Proposition 4.3. Let ^ he a discrete measure ji = ^TTj2xeA^>'^(^ ~ A) such 

that \ax\ < Oo < I and the discrete set A, the support of ji, satisfies the uniform 
discreteness condition (2.5). Let each point A G A move to a new position A' so 
that the mapping $ : A 1— > A' transforms A to another discrete set A', also satisfying 
(2.5), moreover, for any A' e A', 



|a'(A' 



Ae4>-i(A') 



< 



Suppose finally that $ satisfies (4.1). We define the measure n' = 2-k X^a'gA' '^'x^i^— 
A') and the function 

*(^) = i E E «A [In |A' - ,)| - In |A - ^1 + 5?(^(A'-i - A'^))] . (4.13) 
A'eA' Ae4>-i(V) 
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Define also a{z) as aA(0) where X{z) is the point in A , closest to z, and, similarly 
a'{z) = a'{X'{z)) where X'{z) is the point in A' , closest to z. Then the function ^ 
satisfies the Poisson equation A^P = ii' — a,nd 

ll'(z) - a'(z) In |z - A'(z)| + a{z) In \z - \{z)\\ < C\z\^+^ . (4.14) 

The discrete version of Theorem 4.2 is now formulated as following. 

Theorem 4.4. Suppose that /Xq is a measure satisfying conditions of Theorem 4-2- 
Let n be a discrete measure and $ a mapping such that the conditions of Propo- 
sition 4-3 are satisfied. Suppose finally that the union of the sets Ao,A, A' also 
satisfies (2.5). 

Then the measure /xq + — M possesses the potential U{z) subject to 

U{z)-ao{z)ln\z-Xo{z)\+a'{z)lii\z-X'{z)\-a{z)ln\z-X{z)\ > -f'\z\^, < 7' < 7 

(4.15) 

for \z\ large enough, and thus the Pauli operators P_max) resp., P_max have infin- 
itely many zero modes. 

We will give later some examples how Theorems 4.2 and 4.4 can be applied in 
interesting concrete situations. 

Further on wc show that the quadratic growth of the potential and the property 
of the null subspace to be infinite-dimensional are stable also under certain addi- 
tive perturbations of the magnetic field. The perturbation of the field is a signed 
measure /i such that its discrete part /idisc with weights (3\ satisfies conditions of 
Proposition 4.3 and the continuous part ^cont satisfies 

j \\n\w - z\\d\^i^\{w) < B (4.16) 

\w — z\<r^ 

for some B, for all 2; e C (r° is the constant in (2.5)). To evaluate the size of the 
measure yu, the following characteristics will be used: 

w(r) =a;(r,M) = |M|(D(0,r)), 
where, recall, D(0, r) is the disk with center at the origin and radius r, and 



M(r) = / w '^dii{w) 

J ri <\w\<r 



I r\<\w\<r 

for some fixed r\. To such measure /i we associate the potential ^{z) defined as 



1 

2^ 



j \u\l-zlw\dv{w)+ j (^ln|l-i-|+5R(i- + i(i-)'))d.(-) 



\w\<R \w\>R 

(4.17) 

The value of R will be chosen later. Note, moreover, that changing R, we add 

a harmonic function, the real part of a second degree polynomial to ^ ^. For a 
given point z G C, there may be only one point A G A in the r°/2-neighborhood 
of z. If such point exists we denote it by X{z) and set (i{z) = Px{z) and L{z) = 
P{z) In \z — A(z)|, otherwise L{z) is set to be zero. 
Now we formulate our main estimate. 

Proposition 4.5. Under the above conditions the following estimates hold. 
(1) Suppose that for r large enough 

w(r) < C(r2-^),T > 0. (4.18) 

Then 

\^,{z)-Liz)\<C'{\z\ + ir-^ln\z\. 
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(2) Suppose that for r large enough, 

uj{r) < er^ 



\M{r)\ <e,r> Rq. 



Then 



\^,{z)-L{z)\<CeHl/e){\z\ + lf. 
Note that in both cases the estimate can be written as 



(4.19) 

(4.20) 
(4.21) 



|*^(^) - L{z)\ < CQ{\z\)HV^)\zr), 

where Lo{r) is the majorant for Lu{r) in the proposition. 

Before starting to prove the Proposition, we explain, that, in the first case, the 
measure /i is supposed to be rather sparse, so that the |/i|-measure of the disk 
D(0,r) grows slower than the area of the disk. In the second case, the measure of 
the disk may grow as fast as its area, with a small factor, however an additional 
condition (4.20) is imposed. This latter condition is satisfied if the measure /x is 
more or less uniformly distributed in all directions or at least, has zero second 
circular harmonic. On the other hand, if /i is the Lebesgue measure restricted 
to some angle in C not coinciding with the whole plane or a half-plane, then the 
condition (4.20) is violated. 

Proof. We prove the parts (1) and (2) of the proposition simultaneously. Similar 
to the proof of Proposition 4.3, we divide the plane with the disk < R removed 
into the same three parts, = {R < \w\ < ^\z\}, ~ {^|^| < |w| < 2|z|}, and 
^^2 = > ^\z\}- Correspondingly, the integral in (4.17) splits into the sum of 
three integral which we denote by /q, /i, /2- 
For Iq we have 



U{{z/wf)diJ.{w) 



< 



In |1 - z/wldfj^iw) + J ^{z/w)diJ,{w) 
y I In |1 - zM\d\,,\iw) + 1 1 \z/w\d\,,\iw) + \ 




d^i{w) 



.(4.22) 



For w G Oq, we have |2; — > \z\/2 and |1 < |1 — z/'w\ < \2z/w\, and therefore 



1^1/2 

[ \ln\l - -\\dfx(w) < I \n\ — \diiiw) = [ Inl — |da;( 
J w J w J w 



r) 



R 

\z\/2 



= 27r(ln4a;(M)_ln(^)a;(r)) + 27r j ^^^^{r)dr. (4.23) 

R 

The conditions (4.18), resp., (4.19), imply then that the expression in (4.23) is 
majorated by CCj{\z\)). 

The second term on the right-hand side in (4.22) is estimated simply by 

j ^{^)diJL{w) < C\z\ j r-^duj{r) < Cuj{\z\). 

\w\<\\z\ R 
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It is in the last term that the treatment of two cases is different. In the case (1), 

V Oo / ^ 

In the case (2) the estimate by absolute value is not sufficient, so we act differently: 



\ no 



w 



< z 



/ 

mo 



w 



^diJ,{w) 



<\z\'\M{\z\/2)\<e\z\ 



Taken together, the last three inequalities give |/o| < C''^(kl)- 

To estimate Ii, we split the integration set Qi into the disk D = — zj < r°/4} 
and = fii \ D, correspondingly, 7i = /lo + !{■ If there are no points of A in 
the r''/2 -neighborhood of z, so that the measure /i is continuous in D, then Iiq 
is estimated by a constant by (4.16). If there is a (unique) point X{z) £ A in the 
above disk, then this contribution equals p{z) In \z — X{z) \ plus some bounded term 
coming from the continuous part of the measure. In the notation of the Proposition, 
in any case, |/io — L{z)\ < C. 

In the set fl'i, outside the disk D, we again estimate the absolute value of the 
integral by 



|/;| < / I In |1 - z/w\dM{w) + f \z/w\dMw) + 
Jn' Jn' 



z^ /uP'd^{w) 



(4.24) 



The second and the third terms in (4.24) are estimated by l-zp""^ in the same way 
as the similar terms in /q- 

More trouble we have with the first term in (4.24), and the estimates we get for 
this term are the worst ones. We have | In |1 — z/w|| < C + | In | | and the first 
term here contributes with 0(w(|z|)) to I[. We also have 



/ 



zeni,|z-ui|>|z|/2 
So it remains to evaluate 



In 



w 



<Cu;{\z\). 



hi= j \\n\{z-w)/z\\d\n\{w). 

r°/i<\w-z\<\z\/2 

We split the region r°/4 < |w — 2;| < \z\/2 into annuli Uk = {A;r°/4 < — w| < 
{k + l)r"/4}, k=l,2,... (the last one may be somewhat larger). Denote by Sk the 
I /L( I -measure of the annulus Uk- We have, of course, 

Y,Sk<W^'i)<o^{2\z\). 

At the same time, since each annulus Uk can be covered by no more than 8A;7r disks 
with radius r°/4, the condition (4.16) implies that < 'yk, 7 = SttBq. Now, we 
majorize the integral In by the sum of integrals over annuli Uk and estimate it 
from above in the terms of Sk- 

'n<i:in(ll^)s. (4.25) 

The sequence In (^^^ decreases in k. Therefore, if we replace Si by its largest 
possible value 7 with simultaneous decreasing of the rest of Sk to some new non- 
negative values, keeping the sum the same, the sum in (4.25) can only increase. We 
perform the same operation with ^2, making it equal 27, then with ^3 and so on, 
until all 'small' Sk, k < N have their maximal possible values, Sk = k-y, the next 
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one, Sn is smaller than N^, and the rest are zeros. Since the sum of Sk is not 
greater than uj{2\z\), we have ^N{N + l)/2 < ij{2\z\), so N < (2w(2|2;|)/7)V2. in 
this way we have reduced the task of estimating /n to evaluating the sum 



fe=l ^ 



This sum can be estimated by the integral 

.(2a.(2|z|)/7)i/2 /..IN 
Jr" /4 V t'^O J 

After the substitution t = s4:\z\/r°, the integral transforms to 

K2a;(2|^|)/7)i/VNI 

sln{l/s)ds, 



C^\zf [ 
Jo 



which is estimated directly since s\n{l/s)ds — 0(e^ ln(l/e)), e 0. Thus we 
obtain 

/n<Cc.(|z|)ln(y^/|z|). 

Collecting this inequality with previously found estimates for other contributions 
to 7i, we obtain 

\h{z)-L{z)\<Cw{\z\)ln{,MW)/\z\). 
The term I2 is the easiest one. Since \z/w\ < 1/2 in f22, we have 

\ln\l - z/w\ + U{z/w + {z/wf /2)\ < hz/wf. 



f \w\-^d\n\{w) < C\zf [ r-^uj{r)dr, 
Jn2 J2\z\ 



Therefore 

\I2\<1\Z\' 

>J JQ2 •'2\z\ 

which gives the required estimate in both cases. □ 

The Proposition we have proved leads to the following Theorem. 

Theorem 4.6. Let ^0 be a measure subject to conditions of Theorem ^.2. Suppose 
that the perturbation /i satisfies conditions of Proposition 4-5- We suppose also that 
the number e in the case (2) of the Theorem is small enough, so that C ^^elnil / e) < 7, 
where is the constant in (4.21). Suppose finally that the union of the discrete 
sets Aq and A satisfy the condition (2.5), probably, with different r°. Then the 
measure fj,* = iiq + fi possesses a potential 5"* satisfying 

**(z) - a*{z) ln(rf*(z)) > 7*|z|2, 7* > (4.26) 

where d* {z) is the distance from z to the nearest point of the support of the discrete 
part of measure ji* and a* {z) is the measure /z of this point ( or any of such points, 
if there are several of them). In this situation the Pauli operators V-max resp. 
V-EV have infinitely many zero modes 

Although Theorem 3.11 provides us with very general conditions for infiniteness 
of zero modes, it applies only in such cases when there are no arbitrarily large regions 
in the plane where the field is negative - this would destroy the subharmonicity of the 
averaged potential. The perturbation results enable us to establish the infiniteness 
of zero modes in certain situations when such regions are present: they arise as a 
result of perturbations allowed by the theorems in this section. 

The starting point has to be a magnetic field where an exponential estimate 
of the form (4.2) is already known. After this, we can apply allowed types of 
perturbation again and again, as long as at each step the perturbation satisfies the 
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general conditions of this Section, in other words, is weak enough, in the proper 
sense. 

The examples below do not exhaust all kinds of perturbations given by our 
theorems but rather illustrate their possibilities. 

Example 4.7. Let /xq be a measure for which the conditions of Theorem 4.6 are 
satisfied. Let bo a set in C such that the |^o|(f^ H D(0, R)) — 0{R'^~'^) for some 
r > and hq, be the measure restricted to fi. Then for the magnetic field 
^0 ^ BfiQ,, for any positive S the operator Pauli has infinitely many zero modes. 
In fact. Theorem 4.6, part (1) applies here. If _B > 1, the measure /ig — S/io has 
in n the sign opposite to the one of fiQ. The domain ^ can be rather large. In 
particular, for the examples of /xq discussed above, the set may be the domain 
{\X2\ < C(l + Ixil)"^}, with some r < 1, thus admitting rather large negative-field 
regions. Another case is O being the union of disks Dfe of radii Rk tending to 
infinity as A; — > oo, such that the combined area of the disks fitting into the disk 
D(0,i?) is majorated by i?^""^. 

If the field to be perturbed is more regular, even stronger perturbations are 
allowed. 

Example 4.8. Let hq be a regular lattice of Aharonov-Bohm potentials, possibly 
on the background of the constant magnetic field, or a periodic (quasi-periodic) 
measure satisfying Condition 3.6 with nonzero frfcT' '^^^''^^ potential with 
quadratic growth exists. Let fl = 0_ U $1+ be the double angle in the plane, 

n = {argz e (6li,6'2)} U{argz e {9i + 11,62 + n)}, 

and be the restriction of /xq to the angles i'l±. Then, in the notations of Exam- 
ple 4.7, Theorem 4.6, part (2), establishes the required estimate for the potential 
of the measure fi — B/io^ -I- BjiQ and therefore the infiniteness of zero modes for 
the Pauli operator, as soon as the size of the angle, 62 — 6*1, is small enough. This 
smallness guarantees the fulfillment of the condition (4.19), while the symmetry of 
the domain fl leads to fulfillment of (4.20). Instead of symmetry, we can require 
that the domain f2 consists of four angles and is invariant with respect to rotation 
by 7r/2. Then, if the size of the angles is small enough, the perturbation —Bjici 
again satisfies conditions of the perturbation theorems. Of course, one can change 
the field not necessarily in two (four) angles but in any domain symmetric with 
respect to the rotation by tt, as long as the area of the portion of the domain in the 
disk D(0, R) grows not faster than eR^, with e small enough. 

The most efficient pattern for applying the re-arrangement Theorems 4.2 and 
4.4 is the following. Suppose that we have a decomposition of the plane C into the 
union of disjoint sets Q.j so that the diameter of Q.j is not greater than C\zj\'^, r < 1, 
where zj is some point in Qj. Let fi a perturbing measure satisfying conditions of 
Theorems 4.2 or 4.4, such that = 0. Then we can define the mapping as 

$(f2j) = {zj}, so that the whole set flj is mapped into one point. This mapping 
induces the zero measure /i*. The perturbation theorems immediately produce the 
estimates for the potential of the measure fj, together with infiniteness of zero modes. 
This reasoning illustrates that positive and negative parts of the perturbation can 
cancel each other even if they lie rather far apart. 

The sets flj can be chosen rather arbitrarily. The following construction may be 
useful. 

Lemma 4.9. Let t G (0, 1). Then there exists a family of squares Q^, with centers 
at Zj and sides dj covering the plane and having no common interior points, such 
that 

clzj^ <dj<C\zjr, (4.27) 



28 



ROZENBLUM AND SHIROKOV 



with som,e constants c, C , for all squares, except the one containing the origin. 

Proof. We start with the unit square Qo centered at tlie origin, surround Qo by 
eight equal squares which will be called the first layer, so the first nine squares 
form a square with side 3, and condition (4.27) is satisfied, with, say, constants 
c = 1/4, C = 4. Further, inductively, having already a square in the plane filled 
with squares Q^-, so that the inequalities (4.27), we surround this square by squares 
of the same size as in the last layer, if for the new layer (4.27) still holds. In 
the opposite case, if there are N squares in the last layer, the new layer will be 
composed by the squares with the side twice as large as in the last layer (if N is 
even) or with the side 2N/ {N + 1) times larger (if N is odd). One can easily check 
that this construction preserves the inequality (4.27), and repeating it we get the 
required covering. □ 

To show how the above construction works, we consider the case of a continuous 
measure. 

Example 4.10. Let be a continuous measure satisfying the conditions of The- 
orem 4.2 with the following property: for any square Q with center at z, as in 
Lemma 4.9, /i(Q) — -B|Q| < C|zp'^~'^, for some e > with a positive constant 
B, |Q| denoting the area of the square. Wc set ^uq being B times the Lebesgue 
measure; it possesses the required quadratically growing potential and will serve 
as the unperturbed field. The difference, fi — /io will be represented as a sum 
of two terms, /ii + ^2- We set fxi = ^ — iJ,{Qj)/\Qj\dx on the square Q^-, and 
M2 = (m(Qj )/IQj I — B)dx. Then the perturbation satisfies conditions of Theo- 
rem 4.2, with <I> mapping the whole square Qj to the single point Zj G Qj, and /i2 
satisfies conditions of the Theorem 4.6. 

Note that in the situation of the example, the regions where the field points in 
the 'wrong direction', can be very large. In fact, one can choose /i so that for each 
square Qj in our covering, > in the narrow strip near the boundary, with width 
a|2j |^, with an arbitrarily small fixed a, and fi is negative in the main part of the 
square. 

A similar property holds for discrete measures 

Example 4.11. Let A be a regular lattice and ax, A G A is the collection of intensi- 
ties, a\ G [—1/2, 1/2). We suppose that for any square Q with center at z (or only 
for squares Qj constructed in Lemma 4.9), 



for some e > 0. Then, similarly to the previous example, the system of AB solenoids 
with intensities ax placed at the points of A can be obtained by a re-arrangement 
perturbation and an additive perturbation of the regular AB lattice with equal 
intensities. This proves infiniteness of zero modes for the operator V-ev for such 
configuration of the field. Moreover, we can even suppose that initially A is not a 
regular lattice but just a discrete set, in proper sense, almost uniformly distributed 
in the plane. This situation is taken care of by additional perturbation consisting 
in moving the points of A to the points of a regular lattice and then making one 
more additive perturbation to dispose of the points which cannot be moved. We 
do not go into details here. 



^ ax-B\Q\ <C\z 



2t-€ 



B>0 
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